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Abstract: We study neutrinoless double beta decay in left-right symmetric extension of 
the standard model with type I and type II seesaw origin of neutrino masses. Due to the 
enhanced gauge symmetry as well as extended scalar sector, there are several new physics 
sources of neutrinoless double beta decay in this model. Ignoring the left-right gauge boson 
mixing and heavy-light neutrino mixing, we first compute the contributions to neutrinoless 
double beta decay for type I and type II dominant seesaw separately and compare with 
the standard light neutrino contributions. We then repeat the exercise by considering the 
presence of both type I and type II seesaw, having non-negligible contributions to light 
neutrino masses and show the difference in results from individual seesaw cases. Assuming 
the new gauge bosons and scalars to be around a TeV, we constrain different parameters of 
the model including both heavy and light neutrino masses from the requirement of keeping 
the new physics contribution to neutrinoless double beta decay amplitude below the upper 
limit set by the GERDA experiment and also satisfying bounds from lepton flavor violation, 
cosmology and colliders. 
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1 Introduction 

The Standard Model (SM) of particle physics is still the most snccessfnl theory of elementary 
particles and their interactions except gravity. After the discovery of its last missing piece, 
the Higgs boson at the large hadron collider (LHC) in 2012, no convincing sign of new 
physics has appeared nntil the 8 TeV rnn of LHC. Thongh these nnll resnlts are adding 
more feathers to the SM cap, the particle physics commnnity is becoming desperate to 
discover some new physics beyond the SM (BSM). This is dne to the severe inadeqnacies in 
the SM as it fails to explain many experimentally observed phenomena and address some 
theoretical qnestions. One snch observed phenomena is the tiny bnt non-zero nentrino 
masses and large nentrino mixing [1]. Dne to the absence of right handed nentrinos, one 
can not write down a Dirac mass term for the nentrinos whereas the Majorana mass term for 
the nentrinos are disallowed by the gange strnctnre of SM. This keeps the nentrinos massless 
in SM with zero mixing among them which is ruled ont by the recent nentrino experiments 
T2K [2], Donble ChooZ [3], Daya-Bay [4] and RENO [5] . These recent experiments have 
not only made the earlier measnrements of nentrino parameters more precise bnt also led to 
the discovery of non-zero reactor mixing angle @13 which was considered to be (very close) 
zero earlier. The 3cr global fit valnes of nentrino oscillation parameters that have appeared 
in the recent analysis of [ 6 ] and [7] are shown in table 1 and 2 respectively. 
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Parameters 

Normal Hierarchy (NH) 

Inverted Hierarchy (IH) 

Am^i 

I(F5^ 

|Am|i| 

10-3eV2 

sin^ 012 
sin^ 023 

sin^ 013 

<5 

7.02 - 8.09 

2.317-2.607 

0.270 - 0.344 

0.382 - 0.643 

0.0186- 0.0250 

0 - 27r 

7.02 - 8.09 

2.307- 2.590 

0.270 - 0.344 

0.389 - 0.644 

0.0188- 0.0251 

0-27r 


Table 1. Global fit Scr values of neutrino oscillation parameters [6] 


Parameters 

Normal Hierarchy (NH) 

Inverted Hierarchy (IH) 

Am^i 

10-5eV^ 

|Am|i| 

10-3eV^ 

sin^ 012 
sin^ 023 
sin^ 013 
d 

7.11 - 8.18 

2.30-2.65 

0.278 - 0.375 

0.393 - 0.643 

0.0190 - 0.0262 

0-27r 

7.11 - 8.18 

2.20-2.54 

0.278 - 0.375 

0.403 - 0.640 

0.0193- 0.0265 

0-27r 


Table 2. Global fit Scr values of neutrino oscillation parameters [7] 


Although the 3cr range for the leptonic Dirac CP phase <5 is 0 — 27r, there are two possible 
best fit values of it found in the literature: 306° (NH), 254° (IH) [6] and 254° (NH), 266° 
(IH) [7]. There has also been a hint of this Dirac phase to be —7r/2 as reported by [8] 
recently. Although the absolute mass scale of the neutrinos are not yet known, we have an 
upper bound on the sum of absolute neutrino masses from cosmology, given by the Planck 
experiment < 0.23 eV [9]. 

This observation of non-zero but tiny neutrino masses and mixing have led to a sig¬ 
nificant number of research activities in the last few decades in the form of several well 
motivated BSM frameworks. The fact that the neutrino masses are found to lie at least 
twelve order of magnitude lower than the electroweak scale, and the pattern of neutrino 
mixing with large mixing angles is very different from quark mixing with small mixing 
angles, gives the hint that their origin must be different from the physics at electroweak 
scale. The most popular BSM framework explaining the tiny sub-eV neutrino masses is 
the seesaw mechanism which broadly fall into three categories namely, type I [10], type II 
[11, 12] and type III [13], all of which involve the introduction of additional heavy fermion 
or scalar particles into the SM. In generic seesaw models, there exists a hierarchy between 
the electroweak scale and the scale of heavy fermions or scalars required to arrive at the 
suppression for neutrino masses. Therefore in typical seesaw models with order one dimen¬ 
sionless couplings, the additional massive fermion or scalar fields lie at a scale much beyond 
the energies accessible to present experiments like LHC. Bringing these additional particles 
to the TeV ballpark involves the fine-tuning of Yukawa couplings so as to keep the neutrino 
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masses at sub-eV scale. In TeV scale type I seesaw, the Dirac Yukawa couplings have to 
be fine tuned to at least 10“® — 10“^, thereby reducing their production cross sections at 
colliders. TeV scale type II and type III seesaw have slightly better scope of having collider 
signatures due to the presence of electroweak gauge interactions of the additional particles. 
Although the ongoing LHC experiment is still hunting for such new physics signatures, it 
is equally important to look for some independent probe of these seesaw models. One such 
promising arena is the neutrinoless double beta decay (NDBD). For a review, please see 
[14]. 

Neutrinoless double beta decay is a process where a nucleus emits two electrons thereby 
changing its atomic number by two units 

(A, —> (A, Z -\- 2') -\- 2e 

with no neutrinos in the final state. Such a process violates lepton number by two units 
and hence is a probe of Majorana neutrinos, which are predicted by generic seesaw models 
of neutrino masses. There have been serious experimental efforts to detect such a process 
in the last few years. The latest experiments that have improved the lower bound on 
the half-life of NDBD process are KamLAND-Zen [15] and CERDA [16] using Xenon- 
136 and Germanium-76 nuclei respectively. The light SM neutrino contribution to the 
half-life of NDBD can be as small as the lower bound set by these experiments only for 
quasi-degenerate type light neutrino spectrum. Such a mass spectrum of light neutrinos is 
however, tightly constrained from the Planck upper bound on the sum of absolute neutrino 
masses mentioned above. A future observation of NDBD with a half-life close to the present 
lower bound could therefore be a sign of new physics as the Planck upper bound on 
may not allow light SM neutrino explanation for the same. Here we consider a combination 
of type I and type H seesaw mechanisms as the new physics part which, apart from giving 
rise to light neutrino masses, can also give additional contributions to NDBD half-life. We 
also extend the SM gauge symmetry to that of Left-Right Symmetric Models (LRSM) [17] 
where such type I+H seesaw arises naturally. Several earlier works [11, 18] have calculated 
the new physics contributions to NDBD within this model. More recently, the authors of 
[19, 20] studied the new physics contributions to NDBD process for TeV scale LRSM with 
dominant type H seesaw. There have also been several works [21] where type I seesaw limit 
was also included into the computation of NDBD in LRSM. Some more detailed analysis 
incorporating left-right gauge mixing were discussed in the works [22], both in minimal as 
well as non-minimal versions of LRSM. 

In the present work, we consider a scenario where both type I and type H seesaw 
terms can be equally dominating which to our knowledge, was not discussed previously in 
the context of NDBD. Such a combination of type I and type H seesaw together can also 
explain non-zero reactor mixing angle 0i3 as well as CP phases, as discussed in the works 
[23]. Making use of the presence of two equally dominating seesaw terms in the neutrino 
mass formula, we consider the possibility where either type I or type H or both type I 
and type H mass matrices can be arbitrary while structural cancellation [24] between them 
can give rise to the correct light neutrino mass matrix. Instead of considering completely 
arbitrary type I and type H mass matrices, we consider a specific type I seesaw mass matrix: 
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tri-bimaximal (TBM) type which gives @13 = 0 [23]. The type II seesaw mass matrix is then 
constructed in such a way that the correct light neutrino mass matrix is obtained. Assuming 
the left-right symmetry breaking scale to be within the TeV range, we then constrain the 
parameters of the model as well as the relative contribution of individual seesaw terms to 
light neutrino mass formula from the requirement of keeping the NDBD amplitude below 
the upper bound set by experiments. We hnd that the new physics contribution can be 
quite close to or even above this upper bound for some region of parameter space even if 
the light neutrino spectrum is not quasi-degenerate type. This not only allow us to rule 
out some region of parameter space but also increases the possibility of detecting such a 
process in successive run of double beta decay experiments like GERDA. 

This paper is organized as follows. In section 2, we briefly discuss the left-right sym¬ 
metric model and then summarize the origin of neutrino masses in this model in section 3. 
In section 4, we briefly point out the possible new physics sources to neutrinoless double 
beta decay amplitude and discuss them in the limit of type I seesaw dominance, type II 
seesaw dominance and equal dominance of both type I and type II seesaw. In section 5 
we briefly discuss the existing experimental bounds on the masses of new particles in the 
model. In section 6 , we discuss our numerical analysis and hnally conclude in 7. 


2 Left-Right Symmetric Model 

Left-Right Symmetric Model [17] is one of the most popular BSM framework where the 
gauge symmetry of the SM is extended to SU{3)c x SU{2)l x SU{2)r x U{1)b-l- The 


right handed fermions which are singlets under the SU{2)l of SM, transform as doublets 


under SU{2)b, making the presence of right handed neutrinos natural in this model. The 


Higgs doublet of the SM is replaced by a Higgs bidoublet to allow couplings between left and 


right handed fermions, both of which are doublets under SU{2)l and SU{2)b, respectively. 
The enhanced gauge symmetry of the model SU {2)b x 17(1) b-l is broken down to the U{1)y 
of SM by the vacuum expectation value (vev) of additional Higgs scalar, transforming as 
triplet under SU {2)b and having non-zero U{1)b-l charge. This triplet also gives Majorana 
masses to right handed neutrinos, responsible for type I seesaw. The left handed Higgs 
triplet on the other hand, can give tiny Majorana masses to the SM neutrinos through type 
H seesaw mechanism. 

The fermion content of the minimal LRSM is 




Similarly, the Higgs content of the minimal LRSM is 
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r\j 


( 1 , 2 , 2 , 0 ) 


rsj 


(1,1, 3, 2) 


Al 


5l -5t/V‘^) 


~ (1,3,1,2), 


^'r 


where the numbers in brackets correspond to the quantum numbers with respect to the 
gauge group SU{3)c x SU{2)l x SU{2)r x U{1)b-l- hr the symmetry breaking pattern, 
the neutral component of the Higgs triplet A/j acquires a vev to break the gauge symmetry 
of the LRSM into that of the SM and then to the 17(1) of electromagnetism by the vev of 
the neutral component of Higgs bidoublet <1>: 


SU{2)Ly<SU{2)RXU{l)B-L SU{2 )lXU{1)y ^ U{l)em 

The symmetry breaking of SU{2)b x U{1)b-l into the U{1)y of standard model can also 
be achieved at two stages by choosing a non-minimal scalar sector [26]. 


3 Neutrino Mass in LRSM 

The relevant Yukawa couplings which lead to small non-zero neutrino mass are given by 

= Uij^iL^^jR + y'ijhl^^jR + h.c. 

T fij {^Tr ^ i^ 2 ^R(^jR + {R t-)- L)) -|- h.c. (3-1) 

where $ = T 2 ^*T 2 . In the above Yukawa Lagrangian, the indices i,j = 1, 2, 3 correspond to 
the three generations of fermions. The Major ana Yukawa couplings / is same for both left 
and right handed neutrinos because of left-right symmetry (/r = /r). These couplings / 
give rise to the Majorana mass terms of both left handed and right handed neutrinos after 
the triplet Higgs fields Ar^r acquire non-zero vev. These mass terms appear in the seesaw 
formula as discussed below. The resulting seesaw formula in this minimal model can be 
written as 

+ Ml (3.2) 

where the usual type I seesaw term Ml is given by the expression. 

Ml = -ruRRM^l^mlR. (3.3) 

Here turr = yvi y'v 2 is the Dirac neutrino mass matrix, with vi ^2 are the vev’s of the 
neutral components of the Higgs bidoublet. It should be noted that in the framework of 
LRSM, Mrr arises naturally as a result of parity breaking at high energy and both the type 
I and type H terms can be written in terms of Mrr. In LRSM with Higgs triplets, Mrr 
can be expressed as Mrr = vrJr. The first term in equation (3.2) is due to the vev of 
SU{2)r Higgs triplet. Thus, = /rVr and Mrr = JrVr, where vr^r denote the vev’s 
and /l,_r are symmetric 3x3 matrices. The left-right symmetry demands /r = /l = / as 
mentioned above. The induced vev for the left-handed triplet vr can be shown for generic 
LRSM to be 


VR 
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with M\y ~ 80.4 GeV being the weak boson mass such that 


\vl\ « Mw « It's! 

In general 7 is a function of various couplings in the scalar potential of generic LRSM. 
Using the results from Deshpande et ah, (fifth reference in [17]), 7 is given by 

^ I32vI + (iiviV2 + l^zvl .0 

^ {2pi - p3){vl + v^) 

where /3, p are dimensionless parameters of the scalar potential. Without any fine tuning 
7 is expected to be of the order unity (7 ~ 1). However, for TeV scale type I+II seesaw, 7 
has to be fine-tuned as we discuss later. The type II seesaw formula in equation (3.2) can 
now be expressed as 

My = j{Mw/vRfMRR - mLRM~]pnfil^ (3.5) 





Figure 1. Feynman diagrams for Neutrinoless double beta decay due to v—Wl,N — Wr, Ar — Wr 
contributions 


4 NDBD in LRSM 

Due to the presence of several additional vector, scalar and fermionic fields not present in 
the SM, one can have several new physics contributions to the neutrinoless double beta 
decay in LRSM. The corresponding Feynman diagrams are given in earlier works, for ex¬ 
ample [21]. Here we consider only three contributions into account, as shown in figure 1. 
Before discussing these three contributions in detail, we first briefly summarize the different 
possible contributions to NDBD in LRSM below. 

1. The SM contribution comes from the Feynman diagram where the intermediate 
particles are W£ bosons and light neutrinos. The amplitude of this process depends upon 
the leptonic mixing matrix elements and the light neutrino masses. 

2. The light neutrino contribution can come also from the Feynman diagram mediated 
by Wr bosons. The amplitude of this process depends upon the mixing between light and 
heavy neutrinos as well as Wr mass. 
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3. The light neutrino contribution can also come from the Feynman diagram mediated 
by both W£ and W^. The amplitude depends upon the mixing between light hand heavy 
neutrinos, leptonic mixing matrix elements, light neutrino masses and mass. 

4. The heavy right handed neutrino vr contribution can come from the Feynman 
diagrams mediated by WJ^ bosons. The amplitude depends upon the mixing between light 
and heavy neutrinos as well as masses of vr. 

5. The heavy right handed neutrino contribution can also come from the Feynman di¬ 
agrams mediated by bosons. The corresponding amplitude depends upon the elements 
of right handed leptonic mixing matrix and masses of ur. 

6. The heavy right handed neutrino contribution can come from the Feynman diagram 
where the intermediate particles are WJ^ and simultaneously. The amplitude depends 
upon the right handed leptonic mixing elements, mixing between light and heavy neutrinos 
as well as heavy neutrino masses. 

7. The triplet Higgs scalars and can also contribute to neutrinos double beta 
decay through WJ^ and mediation respectively. The amplitude depends upon the 
masses of ^l,r scalars as well as their couplings to leptons. 

For the purpose of this work we consider only two contributions to NDBD in addition 
to the standard light neutrino contribution through exchange mentioned at point 1 
above. These two corresponds to the ones mentioned at point 5 and point 7 above. We 
ignore the Feynman diagrams involving exchange as well as light and heavy 

neutrino mixings. We also ignore the mixing between W£ and bosons, as it is tightly 
constrained from electroweak precision data. Among the triplet exchange diagrams, we only 
consider the one mediated by and A/j as the couplings. The other possible contribution 
through /S.L, W£ is suppressed by the type II seesaw contribution to light neutrino masses 
and hence neglected here. The amplitude of the light neutrino contribution (first Feynman 
diagram in figure 1) considered here is 

PC y^^rniUl^iMujO) (4.1) 

i 

where (following the notation of [25]) is the nuclear matrix element (NME) which 

is a function of the mediating neutrino mass {^) for the NDBD of different nuclei. The 
NME above is written as A4iy{0) as it is independent of neutrino mass for light neutrino 
exchange. In the above expression, Ul is the leptonic mixing matrix which appears in left- 
handed charged current interactions and m* are the masses of light neutrinos for i = 1, 2, 3. 
The contribution from the heavy neutrino and exchange (second Feynman diagram in 
figure 1) can be written as 

Ann oc ^^(4.2) 

where Aiu, Ur is the right handed lepton mixing matrix and Mi are the 

masses of right handed neutrinos for i = 1,2,3. The contribution from W^,A^r exchange 
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(4.3) 


(third Feynman diagram in figure 1) is given by the amplitude 

2 

where — is the additional suppression coming from the scalar propagator with p being 

the average momentum exchange for the process. In the above expression, 

is the (11) element of the right handed neutrino mass matrix. Thus, the standard light 

neutrino contribution can be written as 


NDBD 

ln2 


Gf 


l-M.(0)|^ 


mt 


\ul,m^\^ = GF 


\M.m 


\m 


efr|2 


mt 


(4.4) 


where Gj? is the Fermi coupling constant. Similarly, the new physics (NP) contributions 
considered in this model can be written as 

ln2 m2 M.(0) MiJ 

which can further be simplified to get 


F 


NP 

NDBD 

ln2 


^Wr ^Wr ^Ar 


(4.5) 


or in a shorter notation it can be written as 


F 


NP 

NDBD 

ln2 


G 


l-M.( 0 )|^ 


,efr 


mt 


mjv + 


eff |2 


(4.6) 


where m®®, m^, are the effective neutrino masses corresponding to light neutrino 
(fl = f), heavy neutrino {fr = N) and triplet (A/j) contributions respectively to neutri¬ 
noless double beta decay. 

Our goal in this work is to point out the new physics contribution to NDBD when type 
I and type II seesaw both can be equally dominating. This can be very different from the 
type I or type II dominance cases discussed in earlier works, for example [21]. To show this 
difference we have adopted the simplified approach of earlier work [21] where the Wl — VFr 
and F — N mixing were neglected, resulting in the contributions to NDBD mentioned in 
equation (4.5) being dominant. We leave a more general discussion including all possible 
contributions to a subsequent work. 

Depending on the seesaw mechanism at work, these new physics sources can have 
different contributions to the neutrinoless double beta decay. We discuss type I dominance, 
type II dominance and equally dominant type I and type II seesaw mechanism below with 
reference to their contributions to neutrinoless double beta decay. 


4.1 Dominant Type I Seesaw 


For dominant type I seesaw, the first term on the right hand side of equation (3.5) can be 
neglected. The light neutrino mass can then be written as 


= 




- 8 


VR 




(4.7) 




















where Ur is the diagonalizing matrix of Mrr = /rvh. Multiplying both sides of Mi, by 
Ur in the above equation, we get 


U^MMr = - — U^ruLRURf^^U^mljiUR 
vr ^ 

The most general Dirac neutrino mass matrix vulr is usually diagonalized by a bi-unitary 
transformation = UJ^tulrUr. Assuming Ur = Ur gives 

UlhRUn = 

Thus the light and heavy neutrino mass matrices diagonalized by the same unitary matrix 
Ur in this approximation. Also the light neutrino masses are given by 


rui = 


m 


diag\ ' 
LR 


Mi 


(4.8) 


Thus, in the type I seesaw limit, we can replace the Ur matrix in equation (4.5) by the usual 
neutrino mixing matrix Ur- For diagonal charged lepton mass matrix, the neutrino mixing 
matrix Ur is same as the leptonic mixing matrix, which is known from experimental data. 
The masses of W^, Ar can be hxed within a few TeV, allowed by experimental constraints. 
One of the right handed neutrino masses M* can be hxed within a TeV whereas the others 
can be expressed as ratios of light neutrino masses due to the proportionality oc ^ in 
the type I seesaw limit. For simplicity, the hierarchies of Dirac Yukawa couplings between 
different fermion generations can be neglected so that the heavy neutrino mass ratios can 
be written as mass ratios of light neutrinos, following from equation (4.8). 


4.2 Dominant Type II Seesaw 

If the type I seesaw term is negligible, then the light neutrino is given by the hrst term on the 
right hand side of equation (3.5). Since the light neutrino mass matrix Mi, is proportional to 
the heavy neutrino mass matrix Mrr in this case, the same unitary matrix can diagonalize 
both Mjy and Mrr. Also, in this case the light neutrino masses are directly proportional 
to heavy neutrino masses rrii oc Mi due to the proportionality between respective mass 
matrices. This allows us to express two heavy neutrino mass ratios in terms of light neutrino 
mass ratios. Thus, similar to the type I seesaw dominance case, here also we can write the 
new physics contribution to NDBD in terms of leptonic mixing matrix elements, one of the 
heavy neutrino masses, right handed gauge boson and right handed scalar triplet masses. 


4.3 Combination of Type I and Type II Seesaw 

The new physics contribution to NDBD can be different from the above two cases if type 
I and type II seesaw contributions to light neutrino masses are comparable. Some simple 
relations relating different mass matrices involved in the formula for light neutrino masses 
in LRSM given by equation (3.5) were discussed in [27]. One useful parametrization of the 
Dirac neutrino mass matrix in the presence of type I+II seesaw was studied by the authors 
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of [28]. In another work [29], relations between type I and type II seesaw mass matrices 
were derived by considering the Dirac neutrino mass matrix to be known. If the Dirac 
neutrino mass matrix miR is not known, then we can still choose at least on of the type I 
and type II seesaw mass matrices arbitrarily due to the freedom we have in choosing ttilr 
that appears in the type I seesaw term. After choosing one the seesaw mass matrices, the 
other gets completely fixed if the light neutrino mass matrix is completely known. 

Instead of choosing arbitrary type I and type II seesaw mass matrices, it is really ap¬ 
pealing for us to consider one of these seesaw mass matrices to possess fJ. — t symmetry 
or more specifically, Tri-Bimaximal or TBM type mixing. TBM mixing is a good ap¬ 
proximation to observed neutrino mixing at leading order predicting the mixing angles as 
$12 — 35.3°, 023 = 45° and 0i3 = 0. This TBM mixing matrix discussed widely in the 
literature [30] can be written as 


C^TBM 


('fi 7i M 

_ 111 

V6 Vs 72 
\ V6 Vs V 2 / 


(4.9) 


This TBM type mixing can be accommodated within several discrete flavor symmetry 
models [31]. Since our intention in the present work is to do a phenomenological study of 
equally dominant type I and type II seesaw in the context of NDBD, we do not investigate 
the UV complete flavor symmetry framework of this scenario. The required correction to 
TBM type neutrino mass matrix in order to generate non-zero but small reactor mixing 
angle 0 i 3 can originate from other seesaw terms or corrections from charged lepton sector. 
Here we consider a diagonal charged lepton mass matrix such that the leptonic mixing 
matrix is same as the diagonalizing matrix of the light neutrino mass matrix M^. Although 
either type I or type II seesaw mass matrix can give rise to TBM type neutrino mixing, 
here we consider the type I seesaw mass matrix to be of TBM type mixing whereas type II 
seesaw mass matrix gives rise to the deviations from TBM in order to generate non-zero 0i3. 
Since, type II seesaw term is proportional to the right handed neutrino mass matrix Mrr in 
LRSM as shown in equation (3.5), one can construct Mrr as a deviation from TBM form of 
type I seesaw mass matrix. One can perform this exercise the other way round as well that 
is, assuming the type II seesaw to give rise to TBM type neutrino mass matrix whereas type 
I seesaw gives the necessary correction. Although both of these approaches will produce the 
same light neutrino masses and mixing, they will have different implications in neutrinoless 
double beta decay. Here we consider only the former case, that is type I seesaw mass matrix 
of TBM type leaving the other possibility to future works. 

The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix is related to 
the diagonalizing matrices of neutrino and charged lepton mass matrices Ui,, Ui respectively, 
as 

t^PMNS = UjUu (4.10) 
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The PMNS mixing matrix can be parametrized as 


^^PMNS 


' C12C13 S12C13 

-S12C23 - Cl2S23'Sl3e*'^ C12C23 - 5125235136**^ 
\ 512523 — Cl2C235l3e*^ —C12523 “ 5126235136**^ 


5136"**^ 

523613 

C 23 C 13 


Uuaj 


(4.11) 


where Cij = cos 6 ij, Sij = sin 9ij and 5 is the leptonic Dirac CP phase. The diagonal 
matrix Ums.^ = diag(l, e*", e*^^"*"*^^) contains the Majorana CP phases a, 13 which remain 
undetermined at neutrino oscillation experiments. In the diagonal Ui, the leptonic mixing 
matrix is C/pmns = Ui, which for Uu = C/tbm results in vanishing reactor mixing angle 
013 and the leptonic Dirac CP phase 5. Thus, the type I seesaw mass matrix gives rise to 
vanishing 0 i 3 and 6 whereas type II seesaw mass matrix generates non-zero 0 i 3 and non¬ 
trivial value of Dirac CP phase. Since the diagonalizing matrix of Mi, is Cpmns ^^id that 
of type I mass matrix Ml is DtbM) the formula for light neutrino masses in the presence 
of type I and type II seesaw can be written as 


= Ml^ - (4.12) 


where the Majorana phases are incorporated in the type I seesaw term. In principle, the two 
terms on the right hand side of the above equation can have arbitrary strength provided 
the difference between them gives rise to the correct sub-eV scale light neutrino masses. 
The relative strength of type I and type II seesaw terms can be decided by introducing 
a parameter X such that the diagonal type I seesaw mass matrix can be parametrized as 
^/(diag) _ parameter X can take any numerical values, provided the two 

seesaw terms give rise to correct light neutrino mass matrix, after structural cancellation. 
For a particular value of X, one can construct the type II seesaw mass matrix using the 
above equation (4.12). We denote the symmetric type II seesaw mass matrix as 


M. 


II 


^ Til Ti2 Ti3 \ 
Ti 2 T 22 T 23 
\Ti 3 T 23 T 33 j 


(4.13) 


and using equation (4.12), the type II seesaw mass matrix elements can be derived as shown 
in Appendix A. 

For normal hierarchy, the diagonal mass matrix of the light neutrinos can be written 
as = diag(mi, + Xm^i, -|- Am|^) whereas for inverted hierarchy it can 

be written as = diag(y^m 3 -|- Am 23 — Xm^i, -^/m 3 ~+~Am 2 ^, m 3 ). The mass squared 

differences can be taken from the global fit neutrino oscillation data shown in table 1 and 
table 2 shown above, leaving the lightest neutrino mass as free parameter in Thus, 

the type II seesaw mass matrix can be written in terms of five free parameters: the lightest 
neutrino mass, three leptonic CP phases and the seesaw relative strength factor X. The 
right handed neutrino mass matrix can be found from 


l{Mw/ vr^Mur = mI^ 


which was also shown in equation (3.5). Thus, hxing the dimensionless parameter X allows 
us to calculate the type II seesaw mass matrix in terms of neutrino parameters, which 
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can then be used to find the right handed neutrino mass matrix by fixing 7 and 

vr « Mwn/gR- 


Parameters 

Values (NH) 

Values (IH) 

Am2i 

10-5eV2 

|Am|i| 

10-3eV2 

Sin^ 012 
sin^ 023 
sin^ 013 

P 

Mwl 

7.60 

2.48 

0.323 

0.567 

0.0234 

100 MeV 

80.4 GeV 

7.60 

2.38 

0.323 

0.573 

0.024 

100 MeV 

80.4 GeV 


Table 3. Numerical values of several parameters used in the calculation of m®® for NDBD 
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Figure 2. Standard Model light neutrino contribution to effective neutrino mass which appears in 
NDBD 

5 Lepton Flavor Violation and Collider Bounds 

Lepton flavor violation (LFV) in LRSM were studied in details in previous works including 
[32]. Within this model, there are several possible LFV processes like /i —)• 67 ,// —>• 3e. 
Here we consider ^ 3e process mediated by doubly charged bosons in LRSM. It turns 
out that imposing the experimental bound on this process BR(/i —?■ 3e) < 10 x 10“^^ [33] 
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Figure 3. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type I seesaw dominance and = 3.5 TeV. 


is enough to keep other LFV processes within experimental limit. The branching ratio for 
the /X —>■ 3e process induced by doubly charged bosons A^'*', A^^ is given by [32] 

BR{^^ ^ 3e) = -\h,,Kf f ^ ^ j (5.1) 

where the couplings h are given by 



In equation (5.1), ^ are the masses of A^"^ and in equation (5.2), is the right 

handed leptonic mixing matrix. In a previous work [19], the experimental bound on this 
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Figure 4. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type II seesaw dominance and = 3.5 TeV. 

LFV process was incorporated to restrict where -h ^^2 ■ It 

was found that for most of the parameter space, < 0.1 with My/^ = 3.5 

TeV. Assuming = M^, the above bound will become < 0.1/-v/2- 

However, this bound was calculated only with the assumption that = Kl and hence may 
not be applicable in a general case where both type I and type II seesaw terms contribute 
to light neutrino masses. 

Apart from LFV bounds on the ratio M^®®'™®’'/Ma, there exists other experimental 
bounds on the new particles of LRSM. The most stringent bound on the additional charged 
vector boson Wr comes from the K — K mixing: > 2.5 TeV [34]. Direct searches 

at LHC also put similar constraints on the mass of Wr boson. Dijet resonance search 
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Figure 5. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type I seesaw dominance and Mwf^ = 3.5 TeV, Ma^ = 1 TeV. 

by ATLAS puts a bound > 2.45 TeV at 95% CL [35]. This bound can however be 

relaxed to M\y^ > 2 TeV if gn ~ O.Ggi. There are other bounds on coming from 

other searches in LHC experiments, but they are weaker than the dijet resonance bound 
and hence skipped here. Similarly, the doubly charged scalars also face limits from CMS 
and ATLAS experiments at LHC: 

Ma±± > 445 GeV (409 GeV) for CMS (ATLAS) 

These limits have been put by assuming 100% leptonic branching factions [36]. The heavy 
right handed neutrinos with SU{2)^ gauge interactions are also constrained by direct 
searches at LHC. For example, the search for Wr —>■ IN at ATLAS and CMS constrains 
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Figure 6. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type II seesaw dominance and = 3.5 TeV, Ma„ = I TeV. 

the right handed neutrino masses to be around 1 TeV [37]. All these experimental bounds 
are taken into account in our analysis below. 

6 Numerical Analysis 

We first calculate the standard light neutrino contribution to NDBD (4.4) by evaluating 
the corresponding effective neutrino mass 


16 - 














where Ul = C^pmns is given by equation (4.11). This is effectively the (1,1) (or (ee) in 
ffavor basis) element of the light neutrino mass matrix ^s 

= ufj" = mic?2C?3 + m2sf2C?3e^*" + m3S?3e^*^ (6.1) 


Using the best fit values of three mixing angle and two mass squared differences as shown 
in table 3, one can show the variation of as a function of lightest neutrino mass 
^.lightest = mi(NH), 7713 ( 114 ). This is shown in figure 2 where the Majorana CP phases a, (3 
are varied in the range (0 — 27r), allowed by neutrino oscillation data. It is to be noted 
that the best fit neutrino parameters except 5 in table 3 are taken from [7]. A value of 
Dirac CP phase 6 = —7r/2 has been reported recently by experimental data [8], but the 
standard light neutrino contribution to NDBD is independent of 6 as seen from the above 
expression. It can be seen from figure 2 that the light neutrino contribution can saturate 
the CERDA bound [16] only for higher values of lightest neutrino masses, disallowed by 
the Planck data on sum of absolute neutrino masses [9]. Thus, the Planck constraint on 
the parameter space shown in figure 2 is more strict compared to the CERDA data as far 
as light neutrino contribution to NDBD is concerned. However, the CERDA limit is strong 
enough to rule out many new physics contributions as we discuss below. 

For the new physics contribution discussed in this work, the total effective mass is 


’’T^N+Ar 


2^^ , 


( 6 . 2 ) 

(6.3) 


which, for type I and type II seesaw dominance is calculated following the analysis discussed 
in the previous section. These are shown in figure 3,4, 5 and 6 respectively for two different 
values of masses. For all these cases, the Dirac CP phase 6 is varied in the allowed 3 ct 
range (0 — 27r). Also, Mwr = 3.5 TeV and = 0.1/\/2. Since for either type 

I or type II dominating cases we are considering equality of left and right handed mixing 
matrices (as discussed above), the LEV bound = 0.1/\/2 is applicable here. 

It can be seen from figure 3 that for type I seesaw dominance with Ms = 3.5 TeV, the 
new physics contribution to NDBD is sizable only when rTTughtest > 0.003 eV for normal 
hierarchy. For inverted hierarchy however, the new physics contribution to NDBD is about 
ten times suppressed compared to the standard light neutrino contribution shown in figure 
2. Similarly, for type II seesaw dominance with Ms = 3.5 TeV, CERDA upper limit rules out 
^lightest less than approximately 2 x 10“^ eV and 5 x 10“® eV for NH and IH respectively, 
as seen in figure 4. Of course, these limits are for specific values of 
mentioned above. We have also assumed Ms = M^Tr for the cases shown in figure 3, 4. 
Considering a smaller value of Ms say, 1 TeV and hence a smaller j\^heaviest 
physics contribution to NDBD as 777^^_^^ is inversely proportional to their masses, seen 
from equation (6.2). This is shown in the — TTiughtest plane in figure 5 and 6. 

Now let us consider the interesting scenario where both type I and type II seesaw terms 
are equally dominating. In this case, the simple relation {Ul = Upt) between the diagonal¬ 
izing matrices of left and right handed sectors no longer holds, unlike the cases discussed 
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above. Here we have to diagonalize the right handed neutrino mass matrix explicitly in 
order to calculate its eigenvalues Mj as well as Ur. In the above equation (6.2), Ur is the 
diagonalizing matrix of Mrr which, in terms of type II seesaw mass matrix can be written 
as 

Thus, the right handed neutrino masses are inversely proportional to the dimensionless 
parameter 7 . Therefore, the two new physics contributions in the above equation (6.2) 
have different dependence on 7 as can be seen from the equation below. 




.M2 

U rj2 • oMt ULM 


2 '^Rei I 2 

V. 


Wl 


oc'y 


oc- 

7 


(6.5) 


Although dimensionless parameter 7 should be of order one, in TeV scale type I+II seesaw, 
this has to be extremely fine-tuned. It is straightforward to see that for Mj ~ 1 TeV, 
Mwjj ~ 3 TeV, ~ 0.1 eV, the equation (6.4) gives 7 ~ 10“®. Keeping these parameters 
fixed, if we increase 7 , then the type II seesaw contribution to light neutrino mass will exceed 
the desired range <0.1 eV. This increase can however be compensated by increase in the 
strength (parameterized by X) of type I seesaw contribution so that difference between two 
large contributions from type I and type II seesaw can still give sub-eV scale light neutrino 
masses. 




Figure 7. Constraints on dimensionless parameters X and 7 from CERDA, Planck, LEV and 
perturbative bound for =3.5 TeV and Ms = respectively. 

Since the two new physics contributions have opposite dependences on 7 as seen from 
the equation (6.5) above, it is interesting to see the effects of its variation on m®®. We 
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Figure 8. Constraints on dimensionless parameters X and 7 from GERDA, Planck, LEV and 
perturbative bound for Mw^ = 7 TeV and Mg = vr, respectively. 



Figure 9. Constraints on dimensionless parameters X and 7 from GERDA, Planck, LEV and 
perturbative bound for = 10 TeV and Mg = vr, Mf^^^ respectively. 

first vary the parameters ( 7 , X)continuously and calculate the effective neutrino mass for 
the specific parameter values given in table 3. For each chosen value of we consider 

three possible values of Mg namely, and It is to be noted that the physical 

mass of the doubly charged component of Ajj is given by ~ 2 p 2 V% where p 2 is a 

°R 

dimensionless coupling in the scalar potential. Thus, the maximum mass squared of this 
doubly charged boson is 2 p^^^v^ = where Pmax = is the maximum value 

of dimensionless coupling allowed by perturbative unitarity bound. Also, the Wr boson 
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Figure 10. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type I+II seesaw for = 7 TeV and Ms = 

2 2 

mass is given by ~ Br^r _ Thus, given a specific choice of the value of vr can 

be found out from this formula. We also assume equality of gauge couplings ql = g'R in 
accordance with the left-right symmetry. 

We impose the GERDA bound on the new physics contribution to Planck bound 
on the sum of absolute neutrino masses and LEV bound on BR(^ —)• 3e). We explicitly 
calculate BR(/r —>■ 3e) using the expression given in equation (5.1) without any assumptions 
about the right handed lepton mixing matrix and hence about the couplings hij given in 
equation (5.2). We compute these couplings explicitly for a given right handed neutrino 
mass matrix given by equation (6.4). Since we are calculating the right handed neutrino 
masses for a specific combination of ( 7 , X), we also impose the constraint that the heaviest 
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Figure 11. New physics contribution to effective neutrino mass which appears in NDBD for the 
diagrams shown in figure 1 with type I+II seesaw for = 10 TeV and Ms = 

right handed neutrino mass is where = \/4^, the perturbative upper bound 

on dimensionless couplings. After imposing these four bounds, we finally show the allowed 
parameter space in ( 7 , X) plane. Figure 7 shows the allowed parameter space for = 3.5 
TeV with two different values of Ms = vr, It is clearly seen that only a few points 

with inverted hierarchy are allowed by the constraints imposed. Thus the contribution to 
NDBD will be dominated mostly by the standard light neutrino contribution in this case. 
We see that more and more regions of parameter space are allowed if we consider larger 
values of Myy^ namely, 7 TeV and 10 TeV. They are shown in figure 8 and figure 9 for two 
possible values of Mg = vr, like before. It can be seen from these plots that higher 

values of 7 are allowed by the constraints only if the parameter X is taken to be high. This 
is equivalent to increasing the strength of type I seesaw term in the neutrino mass formula 
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Figure 12. Constraints on the right handed neutrino mass in type I+II seesaw scenario for = 

3.5 TeV and Ms = 


(3.5). 

After finding the allowed parameter space in terms of ( 7 , for three different choices 
of M\Yp,, we show the variation of (for new physics contribntions) with lightest nen- 
trino mass for a particnlar combination of {'y,X). For example, we calcnlate for 
7 = = 10 with M\y^ = 7 TeV and Ms = Its variation with mughtest is 

shown in hgnre 10. It can be seen from this hgnre that most of the parameter space for NH 
is rnled ont for miightest < 0.01 eV whereas most of the parameter space with IH are allowed 
by GERDA npper bonnd on Similarly, we also calcnlate m®® with same combination 
of ( 7 , X) bnt for = 10 TeV and Ms = M™^^. This is shown in hgnre 11. As expected, 
the contribntion to rrf^ is snppressed compared to My/^ = 7 TeV case dne to heavier 
particles mediating NDBD. We choose the other parameter valnes from table 3. The right 
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Figure 13. Constraints on the right handed neutrino mass in type I+II seesaw scenario for = 

7 TeV and Ms = 

handed neutrino masses are explicitly calculated by diagonalizing the right handed neutrino 
mass matrix Mp^. After choosing ( 7 ,A 1 ) as well as the parameters given in table 3, the 
right handed neutrino mass matrix Mpp contains four free parameters (mughtest; o,/S)- 
We vary these parameters continuously and calculate rrf^ whose variation with mughtest 
can be seen from figure 10 and 11 . 

Similarly, we can also constrain the right handed neutrino masses, as for each combina¬ 
tion of ( 7 , X), there exists a corresponding Mpm given by equation (6.4). The lightest right 
handed neutrino mass allowed by all these bounds in type I+II seesaw scenario discussed 
in this work for three different values of Wp mass = 3.5, 7,10 TeV and Ms = 

are shown in figure 12, 13 and 14 respectively. It should be noted that, these bounds have 
been found only by constraining the new physics contribution to NDBD from available 
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Figure 14. Constraints on the right handed neutrino mass in type I+II seesaw scenario for = 

10 TeV and Ms = 

experimental as well as perturbative unitarity bounds. 

7 Results and Discussion 

We have studied the consequences of a combination of type I and type II seesaw on the 
amplitude of neutrinoless double beta decay within the framework of a left-right symmetric 
model. Due to the presence of additional gauge bosons and scalars, this model has several 
new sources of neutrinoless double beta decay. We have considered these extra gauge 
bosons and scalars to be near a few TeV to strengthen their contributions to the amplitude 
of NDBD. For simplicity, we have assumed zero mixing between right handed and left 
handed gauge bosons as well as between heavy and light neutrinos. This results in three 
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important contributions to NDBD shown in figure 1. Within this simplified setup, we 
first show the standard light neutrino contribution to the amplitude of NDBD. After this, 
the new physics contributions to NDBD is calculated assuming either type I or type II 
seesaw to contribute fully to the light neutrino masses and compare with the standard light 
neutrino contribution to NDBD. These are shown in figure 2, 3, 4, 5 and 6 respectively for 
different choices of Ms- We have also incorporated the LFV bounds coming from the 

experimental search for ju, —>■ ey, /r —)• 3e in the calculation. 

After confirming the results of several earlier works in this seesaw limit, we then move 
onto the interesting case where both type I and type II seesaw terms can equally contribute 
to light neutrino masses. This allows us to choose the type I and type II seesaw mass 
matrices arbitrarily, provided they sum up to give the correct light neutrino mass matrix. 
Instead of considering such arbitrary mass matrices, we consider a very specific type I 
seesaw mass matrix which gives rise to tri-bimaximal type neutrino mixing with 612 — 
35.3", 023 = 45" and 0 i 3 = 0. The strength of the type I seesaw term is parametrized by 
a dimensionless parameter X. The type II seesaw mass matrix is them evaluated in terms 
of light neutrino mass matrix, constructed using the best fit neutrino data and the TBM 
form of type I seesaw mass matrix. The type II mass matrix therefore, is written in terms 
of five free parameters: the two unknown Majorana CP phases, one Dirac CP phase, the 
lightest neutrino mass and the dimensionless parameter X. In LRSM, the type II seesaw 
mass matrix is also proportional to the right handed neutrino mass matrix, which allows 
us to construct Mrh in terms of the already derived type II seesaw mass matrix upto a 
constant of proportionality. After choosing the other related parameters as given in table 
3, the constant of proportionality between type II seesaw mass matrix and Mrji is 7 as 
seen from equation (6.4). We show from equation (6.4) that for Ml^ of the order of light 
neutrino masses, and lightest right handed neutrino mass around a TeV, the dimensionless 
parameter 7 has to be fine tuned around 10“®. We scan the parameter space 7 ,X from 
the requirement of keeping new physics contribution to rrf^ below the CERDA upper limit 
as well as to satisfy the LFV and perturbative bounds. The allowed regions are shown in 
figure 7, 8 and 9 for three different values of Wr masses. One can see that larger values of 
7 are allowed only when X is also increased. This behavior can be understood by looking 
at the formula for given by equation (6.5). The first term within brackets on the 
right hand side of (6.5) is inversely proportional to right handed neutrino masses. From 
equation (6.4), the right handed neutrino mass is inversely proportional to 7 and directly 
proportional to type II seesaw term. Therefore, the N — contribution to m®® is directly 
proportional to 7 and inversely proportional to type II seesaw term. From the expressions 
of type II seesaw mass matrix given in appendix A, it can be seen that increases with 
increasing X. Therefore, is directly proportional to 7 and inversely proportional to X 
and hence increase in 7 has to be compensated by an increase in X so as to keep below 
the CERDA upper limit. Similar analysis can also be made for the second term within 
brackets on the right hand side of equation (6.5). 

We then consider a pair of benchmark values of the parameters 7 = 10“®, A = 10 
and calculate the new physics contribution to the effective neutrino mass m®® = m®® of 
NDBD. The variations of m®® with lightest neutrino mass are shown in figure 10 and 11 for 
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Mwj^ = 7,10 TeV respectively. This also allows one to discriminate between light neutrino 
mass hierarchies. Since a particular choice of ( 7 , fixes the right handed neutrino mass 
matrix upto the Majorana CP phases, one Dirac CP phase and mughtest) one can also find 
out the allowed values of right handed neutrino masses for the allowed values of ( 7 , X). We 
have shown the allowed values of lightest right handed neutrino mass Mi with m-iightest in 
figure 12, 13 and 14. We note that with respect to the new physics contribution to NDBD, 
the allowed region of parameter space for = 3.5 TeV is very small, implying that 

the standard light neutrino exchange contribution dominates. As Mw^ is increased, the 
LFV and perturbative bounds become weaker and allows more region of parameter space in 
terms of 7 , X which can give rise to sizable new physics contributions to NDBD. It should 
be noted that the purely type I or purely type II seesaw cases discussed earlier had sizable 
new physics contributions to NDBD. However, they can not be considered as limiting cases 
of the general type I+II seesaw discussed here. This is due to the fact that, 7 = 0 will give 
type I seesaw dominance but with a TBM type mixing matrix, ruled out by experimental 
data. 

We also note that the allowed values of the dimensionless parameter 7 given by equation 
(3.4) are very tiny (~ 10“® if only TeV scale type II seesaw dominates. However, in a 
framework with both type I and type H seesaw, 7 can be larger if the strength of type I 
seesaw Vis also increased. Increasing 7 also reduces the right handed neutrino masses as 
clear from equation (6.4) which is also equivalent to increasing type I seesaw term, which is 
inversely proportional to right handed neutrino masses. Such tiny values of 7 are unnatural 
in most LRSM with TeV scale type H seesaw mechanism, and demands the role of some new 
physics behind it. Due to necessity of such unnatural fine-tuning of 7 , there have been many 
studies of LRSM where the terms in the scalar potential leading to the f3 terms appearing 
in the expression for 7 in equation (3.4) are removed by imposing some symmetries, leading 
to Vl = 0 from the minimization of the scalar potential. This will give rise to a purely type 
I seesaw framework with TeV scale left-right symmetry. Such possibilities were discussed 
in the last reference of [17] and also in [38]. 

We have considered a very simplified picture of neutrinoless double beta decay in left- 
right symmetric model ignoring the contributions from left-right mixing as well as heavy- 
light neutrino mixing. However, we have pointed out the differences in for individual 
seesaw dominance and equal dominance of both type I and type H seesaw, and constrained 
the model parameters in a way not considered before. However, a more detailed analysis 
taking into account all the new physics contributions to NDBD in LRSM should be pursued 
to give a more general conclusion. We also did not construct the UV complete flavor 
symmetry framework giving rise to the TBM form of type I seesaw mass matrix. It is 
undoubtedly a non-trivial exercise to implement discrete flavor symmetries in left-right 
symmetric models due to the difference in gauge structure and lepton representations from 
that in the standard model. Very recently such a model building work appeared in [39]. 
We leave such a detailed analysis to a subsequent work. 
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A Elements of Type II Seesaw Mass Matrix 

Til = (^C? 2 C ?3 + mi - ^6^*“ {-X - 3s? 2C?3) + S?3e^*^m3 (A.l) 

Ti 2 = “ 3ci2Ci3Si 2C23 — 3e*'^Si3S23Ci2Cl3)mi — A — 3 si2Ci2Ci3C23 

+3e*^Si2Si3S23Ci3)m2) + ^(3si3S23Ci3e*^^^+'^^m3) (A.2) 

Ti3 = + 3si 2S23C12C13 “ 3e*'^Si3Ci2Ci3C23)mi — A + 3 S 12 S 23 C 12 C 13 

+3e*'^Si2Si3Ci3C23)m2) + ^(3si3Ci3C23e*^^^+'^^m3) (A.3) 


T22 = 


iS A\ 

'S12C23 + e Sl3'S23Cl2j + "^ ) "^1 + 


jA \ ^ a \ 

C12C23 — e Si 2S13'S23j 1 e m2 


(A.4) 


T23 = 


JS, 


JS, 


X 


C23S12 + e Sl3'S23Cl2 j —'S12S23 + e S13C12C23j + "^ ) "^1 

C12'S23 + e*^Sl2'Sl3C23) (ci2C23 “ e*^'Si2Si3S23)-e^*“m2 

+ fs23C?3C23 - e^*(^+'^)m3 


(A.5) 



Sl2'S23 + e*'^Sl3Cl2C23 



mi + 


^(ci 2S23 + e*^Sl2Sl3C23) + y^ e^*“m2 

+ (44 + f) 

(A.6) 
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